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Small systems (of interest in the areas of nanophysics, quantum information, etc.) are particularly 
vulnerable to environmental effects. Thus, we determine various thermodynamic functions for an 
oscillator in an arbitrary heat bath at arbitrary temperatures. Explicit results are presented for 
the most commonly discussed heat bath models: Ohmic, single relaxation time and blackbody 
radiation. 
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1 I. INTRODUCTION 

Heat bath models are of topical interest especially in areas such as quantum information and nanophysics. Thus, 
it is important to develop realistic calculations that can be used to make contact with experiments. Here, we wish to 
examine the effects of a heat bath on various thermodynamic functions such as entropy, partition function, average 
, energy, specific heat and heat capacity. Our starting point is based on an exact result which we have previously derived 
for the free energy of an oscillator in an arbitrary heat bath, in terms of a single integral involving the generalized 
susceptibility [l[ arising from the associated quantum Langevin equation [2] • This result was used in a series of papers 
to obtain free energy shifts of atomic levels in a blackbody radiation field; the effect of a heat bath on the magnetic 
moment of an electron gas [3|, based on a generalization of our previous work to include a magnetic field; a proof 
that the third (Nernst's) law of thermodynamics is valid in the presence of a heat bath 4] and a demonstration that 
a supposed violation of the second law is only apparent [f| . Thus, because of its wide applicability, we are motivated 
to systematically develop explicit results for the most commonly discussed heat bath models. Hence, in Sec. II, 
we review our starting-point Hamiltonian describing an oscillator in an arbitrary heat bath at temperature T which 
enabled us to obtain the equation of motion of the oscillator in terms of a quantum Langevin equation which, in turn 
led us to an exact expression for the free energy of an oscillator in an arbitrary heat bath. Next, we use this general 
result to consider in detail the most commonly discussed heat bath models, obtaining results for the free energy F(T) 
which incorporates the Ohmic, single relaxation time and blackbody radiation models in a form which is very similar 
for all cases, involving the Stieltjes J-function, whose properties we present in Appendix A. This enables us to obtain a 
simple expression, in terms of the J-function for the free energy F(T) which incorporates the Ohmic, single relaxation 
time and blackbody radiation models. This expression for F(T) is then used to obtain explicit results, both for low 
temperature and high temperature, for various thermodynamic functions such as the specific heat, the energy and the 
heat capacity; these results are presented in III (for the Ohmic model) and IV (for the single relaxation and blackbody 
radiation models). Results for the T = case are given in Sec. V. We conclude with a brief summary and discussion 
in Sec. VI. 
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II. FREE ENERGY 



The most general coupling of a quantum particle coupled to a linear passive heat bath is equivalent to an 
independent-oscillator model p], [|J , which is described by the Hamiltonian 



(2.1) 



Here x and p are the particle coordinate and momentum operators and V(x) is the potential energy of an external 
force. The jth independent oscillator has coordinate qj and momentum pj and the generality of the model arises from 
the infinity of oscillators with an arbitrary choice of the mass rrij and frequency ujj for each. 
Use of the Heisenberg equations of motion leads to the quantum Langevin equation 

mx+ f dt'fi{t - t')x(t') + V'(x) = F(t), (2.2) 

J — oo 

where fi(t) is the so-called memory function. F(t) is the random (fluctuation or noise) operator force with mean 
(F(t)) — 0. The quantities n(t) and F(t) describe the properties of the heat bath and are independent of the external 
force. 

In the particular case of an oscillator potential 

V(x) = ^Kx 2 = imc^x 2 . (2.3) 
Substituting ()2.3|) into (|2.2| enables us to obtain the explicit solution 

x{t) = I dt'G(t-t')F(t'), (2.4) 

J — oo 

where G is the Green function. The Green function vanishes for negative times and its Fourier transform, 



a(w) - / dte lut G{t) (2.5) 
J o 

is the familiar response function (generalized susceptibility). This is given by 

<*(z) = 5 1 , (2.6) 

where jl(z) is the Fourier transform of the memory function: 

/>oo 

£(«) = / dtfi(t)e lzt . (2.7) 
Jq 

Note that fi>(z) and, hence, also a(z) are analytic in the upper half plane. 

The system of an oscillator coupled to a heat bath in thermal equilibrium at temperature T has a well-defined free 
energy. The free energy ascribed to the oscillator, F(T), is given by the free energy of the system minus the free 
energy of the heat bath in the absence of the oscillator. This calculation was carried out by two different methods 
d, @| leading to the " remarkable formula" 

fW-lf^l ^*^' }, (2.8, 



7T 







where f(u>,T) is the free energy of a single oscillator of frequency u>, given by 

/(w,T) = kTlog[l-exp(-hu/kT)]. (2.9) 
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Here the zero-point contribution (hui/2) has been omitted, but in a brief Sec. V we remark upon this contribution. 
We have referred to (2.8) as a "remarkable formula" @, Q, in the sense that it displays a non-trivial dependence 
on the temperature T, in contrast with the corresponding classical formula. We have now all the basic tools at our 
disposal and we proceed to consider three cases of interest: 



fl(z) 
AO) 

Here £ is the Ohmic friction constant, while r is the relaxation time. It is generally assumed that the relaxation time 
is small in the sense that r -C C/ m - I n the QED case, fl is a high frequency cutoff characterizing the electron form 
factor (fi — * oo corresponds to a point electron). The susceptibility for all three cases may be combined in a single 
expression 



C, Ohmic, 

c 

; — , Single relaxation time, 

1 — izr 

2e 2 ztt 2 

— — Quantum electrodynamics (QED). (2-10) 

3c d (z + iil) 



in 

-m(z + iQ')(z 2 + ijz — w 2 ) 



For the single relaxation time model 



_ i _ i C__ n' 2 + 1 n' + to 2 k_ 2 n' 

T ~fi~f>' + 7 ' m" 7 (0' + 7 ) 2 ' m~ °SV + Y ' 



The Ohmic model corresponds to the limit of ft' — > oo, in which case r — > 0, £/m — > 7 and if/m — > to 2 . For the QED 
model 



1_ 1 7 K _ 2 Q' M _ (luq + 7fi')(n' + 7) 

SI ^ SI 7 + ^2"' M ~ Wo f2' + 7' ^" ~ ^pT 



where m is the bare mass and 



2e 2 fl 

M = m + — r 2.14 
3CT 

is the renormalized (observed) mass. In this QED case, the limit fl' — > 00 corresponds to the largest value of the 
cutoff £1 consistent with a positive bare mass, that is, in this limit m — 0, K — Mlo 2 and fl = l/r c , where 

' = 6 x 1(T 24 S . (2.15) 



3Mc 3 

With the general form (2.11) the free energy (|2.8|) can be written 



F(T) = — r^log(l- e -^/ fcT )f-^^ + ^— + . „ u;2 + ajg „ - V (2.16) 
v 1 ir J bv y \ u; 2 + Q 2 cj 2 + f7' 2 (uj 2 - uj 2 ) 2 +j 2 lo 2 J y ' 

We use partial fractions in the third term by introducing 



2i=2+ iw i' z l = 2~ iuJl > Wl = Y^o - "J - ) ( 2 - 17 ) 

and we note that, for the overdamped case [(7/2) > wo], &i is imaginary, in which case z\ — 5^ — |u>i| and = ^ + 
Hence 



+- 



+ 



uj z + zf uj z + z 



*2 



I y 2nkT 

where J(z) is the Sticltjes J-function 



-•rlA-A-^)-^)!, 



'2nkT' 



'2irkT' 



(2.18) 



J(z) 



JO 



(itlog(l-e- 27rt ) 



t 2 + z 2 ' 



Imz > 0. 



(2.19) 



In the next two sections, we consider the three specific models separately and in detail. For this purpose, we make 
extensive use of the J-function, whose properties are discussed in detail in Appendix A. 



III. OHMIC MODEL 



Here 



F(T) 



-kT 



'2nkT' 

where in the expression (2.17) for z\ and z\ we put u> — \jKjm and 7 = C/m. 



(3.1) 



A. Low temperature expansion (kT « htoo) 
In the low temperature case we use the asymptotic expansion (A5) for J. With this we obtain for the free energy 



F(T) = -kTj2c n 



The entropy is 



n=0 



7r(fcT) 2 7 7T 3 (fcT) 4 7 (3^ 2 - 7 2 ) 8^ 5 (fcr) 6 7 (5o> 4 - 5 7 2 o> 2 + 7 4 ) 



6fia> 2 



315ft 5 c^° 



(3.2) 



S(T) 



= k 



dF(T) 
dT 

irkTj 4^ 3 (fcr) 3 7 (3co 2 -7 2 ) 16^ 5 (fcT) 5 7 (5w 4 - 5 7 2 w 2 +7 4 ) 



+ 



105h 5 t 
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(3.3) 



The energy is 



+ 



167r 5 (fcT) 6 7 (5^ - 5 7 2 ^ + 7 4 ) 87T 5 (fcT) 6 7 (5w 4 - 5 7 2 w 2 + 7 4 



105ft 5 wJ° 



U{T) = F + TS 



^(fcT) 2 7 , 7r 3 (fcr) 4 7 (3c 2 - 7 2 ) 87r 5 (fcT) 6 7(5a> 4 - 57 2 ^ 2 +7 4 ) , 



(3.4) 
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The specific heat is 

.as 
W 

'irkTj 47r 3 (fcT) 3 7 (3w 2 - 7 2 ) 16^ 5 (fcT) 5 7 (5wg - 5 7 2 ^ 2 + 



(3.5) 



As a check, we note that the leading term in (3.2) agrees with the result obtained by us in p[ while the leading term 
in (3.3) agrees with our earlier results Q as well as a recent result of Hanggi and Ingold Q. In addition, the first two 
terms in (3.5) agree with the results obtained in Q. 

B. High temperature expansion (kT >> hojo) 
In the high temperature case we use the small argument expansion (A4) for J, with the result 

„/„s kT fvy , 2nkT fvjj\ 7 fvy , 

F(T) = -fcTlog— --Ilog- iarccos-i---I 1- 7E 

JUjJq ZlT TOjJq TT ZlOq ZlT 



^ ' n \2nkT 



cos ( n arccos — — ) . (3.6) 
1 2w 7 V ' 



2n 



n=2 

As a check we consider the uncoupled oscillator. Forming the limit 7 — > 0, we find 

n— 1 v / 

= fcTlog(l- e - nMo/kT ), (3.7) 
which is the familiar result (2.9) for the uncoupled oscillator. Here we have used the formula [8] 

log(l - e"«) = logz - \ + E(-)" +1 ^ ■ (3.8) 

71=1 

Returning to the expansion (3.6), we obtain explicit expressions for the first few terms, 

cvm !,ti kT ^7 2nkT Huj 1 7 fvy 

F(T) = -kT log- — log — arccos — (1 - -f E ) 

Thoq Zir JujJq it Aloq Zir 



n 2 (2c 2 - 7 2 ). C(3)^ 3 7(3^~7 2 ) 
48/fcT 24tt 3 (A;T) 2 

With this, the entropy, energy, and specific heat are given, respectively, by 

m - -SP 

, „ kT . fvy 
= fc(log_ + l)' 



o; V\ \n( n_1 X( n ) / ^0 V / 7 \ 

-2fc > — — — — — cos narccos- — 

^ n \2irkT ) \ 2u> J 



n=2 



fc(l g— + 1) 



(3.9) 



hujQ 2ttT 

h^t-^) C(3) /l 3 7(3^-7 2 ) f3101 
48(fcT) 2 12^ 3 (fcT) 3 ' ld Uj 
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U{T) = F + TS 



hr< ( 2nkT \ hwi 7 
kT - — log — j E arccos - — 

2lT \ TVjJq J TT ZUJq 

-2*Tg(-rC(») (^J cos („ arccos 
fcT - — log — 7ij arccos - — 

111 \ fUjjQ J 7T ZUJq 

n 2 (2u;g-7 2 ) C(3) ^ 3 7(3^ 2 - 7 2 ) 
24fcT 8tt 3 (fcT) 2 



(3.11) 



and 



c(r) = t§ 

_ % ,^^,.C(3)^M-j!) m2 i 

" fc 2^T 24(fcT) 2 + 4tt3 (JfeT)3 ' 

Note that all these results apply to the overdamped case with the prescription 

7 7 l^il 
wi arccos- ► log(- ). 

2u>Q 2Ll>q UIq 

Also, we again have a check in that the first three terms in the specific heat agree with the results obtained in Q for 
the Ohmic model. 

IV. SINGLE RELAXATION TIME AND NONRELATIVISTIC QED MODELS 

The free energy is now of the general form (|2.18[) . which can be written 



F(T) — FohmicCO + kT 



J( 2 7 rfcT ) J( 2 7 rfcT ) 



(4.1) 



We argue that f2 and f2' will always be large compared with kT, so it is appropriate to use the low temperature 
expansion and then only the first term. The result is 

F(T) = F 0hmic (T) + «^ ( I - i- V (4.2) 



6h \n n'. 

For the single relaxation time case, it is clear from (2.12) that the second term in (4.2) is very small so that the results 
in this case are essentially the same as for the Ohmic case. However, for the QED case 

so 



tt(A:T) 2 7 

F(T) = Fo^(T) + -^I. (4.4) 



"0 
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A. Low temperature expansion (kT « huio) 

The second term in (4.4) is exactly the negative of the leading term in the low temperature expansion (kT « Hloq) 
for the Ohmic case, given in (3.2). In other words, the T 2 term vanishes and the leading term is the T term. The 
result is that 



-Fqed(T) = - 



7r 3 (fcT) 4 7 (3w 2 - 7 2 ) 87T 5 (fcT) 6 7 (5w 4 - 5 7 2 cj 2 + 7 4 ) 



A5h 3 u 



315^ 5 ( 



,10 



(4.5) 



Sqed(T) = k 



47r 3 (/fcT) 3 7 (3w 2 - 7 2 ) 167r 5 (fcT) 5 7 (5w$ - 5 7 2 w 2 + 7 4 ) 



105^ 5 ( 



,10 



(4.6) 



and 



C/qed(T) 



^ 3 (fcT) 4 7 (3w 2 - 7 2 ) 87r 5 (fcT) 6 7 (5w 4 - 5 7 2 ^ 2 + 7 4 ) 



15ft 3 ^ 



Cqed(T) = k 



47r 3 (fcT) 3 7 (3w 2 - 7 2 ) 167r 5 (fcT) 5 7 (5w 4 - 5 7 2 w 2 + 7 4 ) 



(4.7) 



(4.8) 



We note that, in the large cut-off limit 0], 7 = w§r e , where r e is given in (2.15). In this limit and with 7 << ujo, we 
have a check in that the leading terms in the free energy and the entropy agree with the results obtained earlier by 
us 0. 



B. High Temperature Expansion (kT >> hu>o) 
With the high temperature expansion (3.9) for Fohmic we find from the general expression (4.4) 



, , kT tt(A;T) 2 7 

F QED (T) = -kTlo g — + ^-^ + 



(4.9) 



wr) = ,( (log - +1) -4gi + ...}, 



U qed (T) = kT ■ 



(4.10) 
(4.11) 



and 



, ( 4 7 r 3 (fcT) 3 7 (3c 2 - 7 2 

c Q™( T ) = k [ ^4 

We note that these results agree with the corresponding results in [l[ . 



(4.12) 



8 



V. ZERO-POINT ENERGY 

Since F = U + TS, the zero-point free energy is always identical with the zero-point energy. The zero-point free 
energy is obtained by replacing f(u>,T) — > hu/2 in the formula (|2.8p . The resulting expression diverges for the QED 
model, whatever the cutoff. For the single relaxation time model it is finite for finite relaxation time, 

. . h f°° , fluj Q'co z\uj z\u 

( F > *ero-point = ^ J Q ^^2 + ^2 + ^2 + q,2 + u 2 + z 2 + ^2 + z *2 I 

= — {ft log ft — ft' log ft' — z\ log z\ — z* log } 

Z7T 

ft ft' + 7 ft' + 7 7 

= — {ft' log— — 1 +l\og - + 2w ia rccos-^}. (5.1) 

Zir il ujq 2ujq 

In the Ohmic limit this is logarithmically divergent, 

(-^zero-point ~ ^{tU " log Wor) + 2w x arccos £-}. (5.2) 
VI. CONCLUSIONS 

Motivated by the fact that environmental effects play an important role in many topical areas of physics, where 
dissipation and fluctuation effects often play a significant role, we have presented an exact calculation of quantum 
thermodynamic functions for an oscillator in an arbitrary heat bath at arbitrary temperatures. Explicit results 
were obtained for both high and low temperatures. Since we are dealing with non-additivity of entropies [H, we 
use a method based on (2.8), which is an exact result for the free energy of an oscillator which takes into account 
interaction effects. In the Introduction, we have already given examples of its application 3, 4, 5]. However, there are 
many other possible topics where such results are likely to be applicable. For example, Jordan and Buttiker [TJ] have 
demonstrated the relation between entanglement (due to the heat bath) and energy fluctuations and concluded that 
large entanglement implies large energy fluctuations. Since their work was confined to zero temperature, it would 
be of interest to extend it to non-zero temperatures. In a similar vein, the decrease of the coherence length of an 
Aharonov-Bohm-like interferometer due to interaction with the environment was examined but again it was confined 
to zero temperature fl5j . 

Finally, we turn to a very different area where thermodynamic considerations play a vital role i.e. the study of 
black holes. Following the remarkable results of Bekenstein and Hawking [l6l |. there has been continuing interest in 
developing a microscopic theory for the entropy of a black hole and, in particular, the fact that it depends on the 
area of the event horizon. As an example, we mention the work of Bombclli et al. [171 ] and Srednicki [18| |. where the 
use of partial traces and reduced density matrices played a crucial role. Since, in general, such techniques lead to 
results different from those obtained by the method discussed above, we feel that it would be worthwhile to apply our 
approach to the study of the thermodynamic properties of black holes. 
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APPENDIX A: THE STIELTJES J-FUNCTION 

The Stieltjes J- function is introduced by the integral: fioj 



1 f°° z 
J(z) = — dt\og(l - e- 2 ^)^— 5 , Re{z}>0. (Al) 

7T Jq Z + t 

The imaginary axis is a "natural boundary" of J(z). That is, the analytic continuation of J(z) into the left half plane 
is not given by the analytic continuation of the integral. 
This analytic continuation is based on the identity: fiol] 

J{z) = log[r(z + 1)] - log - (z + i) log(z) + z. (A2) 

Since T(z) is analytic in the entire plane except for poles at z = 0, — 1, —2, • • • , we can use this form throughout the 
z-plane cut along the negative real axis. It is then a simple matter to show that the continued J-function is given by 

J(ze ±M ) = - J{z) - log(l - e^ 2 ™ 2 ), Re(z) > 0. (A3) 
For \z\ < 1, we have the expansion 

J(z) = -logV2^-(z + l/2)log z + z- lE z + Y z n (A4) 

n=1 

where 7^=0.5772157 is Euler's constant and C( n ) is t ne Riemann ^-function. For large z we have the asymptotic 
expansion: [ll| 

J(Z) = S (2n+l)(2n + 2) ~z^ ' (A5) 



where the Bernoulli numbers are 



B2 I' B4 -~^o> B "-h' B8 ~~^o> 

5 R _ 691 7 3617 

66' 2730' 6' 510 ' 

43867 174611 854513 
B x % - B 20 - B 22 . (A6) 



798 ' " 330 ' 138 

Very useful for numerical computation is the Lanczos formula [T3 , Il3j : 



1 4- 4- 1 1 N d 

J(z) = (z + -)log-— --- I- 7 _-+log{d + y;-5-}, ^z>0, (A7) 

2 z 2 ' z + n 

n— 1 



where, for N = 6, 7 = 5, and 



d = 1.000000000190015, di = 76.18009172947146, 
d 2 = -86.50532032941677, d 3 = 24.01409824083091, 
di = -1.231739572450155, d 5 = 0.001208650973866179, 

d 6 = -0.000005395239384953. (A8) 
The numerical error is small (less than a part per billion) everwhere in the right half plane. 
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